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ABSTRACT 



Aims. We study the massive stars in OB associations and their surrounding interstellar medium environment, using a population 
synthesis code. 

Methods. We developed a new population synthesis code for groups of massive stars, where we model the emission of different 
forms of energy and matter from the stars of the association. In particular, the ejection of the two radioactive isotopes 26 Al and 60 Fe is 
followed, as well as the emission of hydrogen ionizing photons, and the kinetic energy of the stellar winds and supernova explosions. 
We investigate various alternative astrophysical inputs and the resulting output sensitivities, especially effects due to the inclusion of 
rotation in stellar models. As the aim of the code is the application to relatively small populations of massive stars, special care is 
taken to address their statistical properties. Our code incorporates both analytical statistical methods applicable to small populations, 
as well as extensive Monte Carlo simulations. 

Results. We find that the inclusion of rotation in the stellar models has a large impact on the interactions between OB associations 
and their surrounding interstellar medium. The emission of 26 Al in the stellar winds is strongly enhanced, compared to non-rotating 
models with the same mass-loss prescription. This compensates the recent reductions in the estimates of mass-loss rates of massive 
stars due to the effects of clumping. Despite the lower mass-loss rates, the power of the winds is actually enhanced for rotating stellar 
models. The supernova power (kinetic energy of their ejecta) is decreased due to longer lifetimes of rotating stars, and therefore the 
wind power dominates over supernova power for the first 6 Myr after a burst of star-formation. For populations typical of nearby 
star-forming regions, the statistical uncertainties are large and clearly non-Gaussian. 

Key words. Stars: abundances, early type, winds, outflows - ISM: abundances - Gamma rays: observations 



1. Introduction 

Massive stars dynamically shape the interstell ar medium 
aroun d them on timescales of a few Myr (s ee e.g. |Lozinskaval 
2001) . Due to their high mass loss rates dCastor et al.l 11975; 
Lamers & Mortonll976tlBarlow et alJl981 | ) and terminal veloc- 
ities of their winds dHow arth & Priniall989tlLamers et al.lll995l) 
and their supernova explosions dJones| [1998) they are important 
sources of mechanical power causing interstellar medium 
(ISM) turbulence and the formation o f shells and cavities (e.g. 
van der Huchll Il987t iLeitherer et ail Il992t iMaeder & Contll 
19941) . Their large UV luminosities, they furthermore cause the 
surrounding medium to be photo ionized to significant distances 
dPanagialll973tlVacca et alJll996l) . 

The mixing of ejecta from young stars into the interstellar 
medium is an important process in the interplay between star 
formation and galaxy evolution. A unique window into these 
processes is provided by the radioactive isotope 26 Al, traced 
by its y-ray decay line at 1808.63 keV. With a mean lifetime 
of ~1 Myr it is a long-term tracer of nucleosynthesis for 



populatio ns of sources able to e ject it sufficiently fast after 
synthesis dPrantzos & Diehllll996l) . The COMPTEL instrument 
aboard the Compton Gamma Ray Observatory (CGRO) has 
mapped the 26 Al emission at 1809 keV in a 9-year full sky 
survey. From a comparison between the image morphology and 
the known spiral arm tangents and regions of star formation, 
one can deduce that massive stars dominate the Galactic 26 Al 
production, and that the co ntribution from no v ae and AGB stars 
must b e relative l y min o r dDiehl et al. 1995 : iKnodlseder et al ' 
1999; Pliischkel 120011: iPrantzos & Diehll (l996t iDiehl et al 
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2006). Localized groups of massive stars in star-forming regions 
such as Cygnus and Orion emit clear 26 Al signals, and these 
regions provide an important contribution to the total amount of 
26 Al present in the Milky Way, and it is therefore important to 
understand these basic building blocks. Typically, massive stars 
eject a few 10~ 5 M o of 2 6 Al through their winds and supernova 
(SN) explosions (see e.g. lLimongi & C hieffi 2006), and the total 
mass of 26 Al in t he Milky Way is estimated to be 2.8±O.8M 
dDiehl et alj|2006l) . 

A similar tracer is the isotope 60 Fe, ob s erved by its 1173 
keV and 1333 keV decay lines dSmithl 12004 lHarris et al.l 
120051) . This isotope is also emitted in the supernova explosions 
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dLimongi & Chieffil 120061) . The lifetime of 60 Fe has recently 
been revised to ~3.8 Myr, from the previous estimate of ~2 
Myr (G. Rugel, private communication). While the lines are 
significantly weaker than the 26 Al line, it is possible to detect 
them when integrating over large regions of the Milky Way. 

We have developed a population synthesis code that follows 
the evolution of massive stars and computes the ejection of 26 Al 
and 60 Fe from a star-forming region. To facilitate studies of the 
dynamics of the surrounding environment, the code also com- 
putes the kinetic energy and the mass ejected from the winds 
and the SN explosi ons, as well as t he ion i zing flux. In contrast 
to St arburst99 ( iLeitherer et al.l 119991: IVazquez & Leithererl 
120051) . we focus on the study of relatively small, nearby star- 
forming regions in the Milky Way. Due to the limited number of 
massive stars in thsee individual regions, stochastic processes 
can cause their behaviour to deviate strongly from the mean 
results obtained for large populations. Because of this, the code 
is designed to not only calculate the average behaviour from 
an analytical integration over the initial mass function (IMF), 
but also the expected distribution of output quantities for finite 
populations. This is both done using analytical approximations 
and Monte Carlo si mulations. To do this we im plemented the 
statistical methods o fC ervifio & Luridianal(l2006l) . These enable 
us to determine when Monte Carlo simulations are mandatory 
in order to properly describe the distributions, and how many 
simulations are needed to obtain statistically robust results. 

Here we describe the population synthesis code, and dis- 
cuss the dependence on the various physical input models. We 
compare our results to the predictions o f Starburst9 9 for ki- 
netic luminosity and UV flux, and to ICervino et al. (2000a) 
for radioactive isotope production. We discuss the simulated 
stellar populations in the context of the inter-stellar medium 
(ISM), with emphasis on the change s induced by recent stel- 
lar evolution tracks for massive stars dLimongi & Chieffill2006t 
iMevnet & Maederll2005l) . 

2. Population synthesis of massive stars 

The synthesis code calculates time profiles of the emission of 
radioactive elements 26 Al and 60 Fe from a population of massive 
stars, for an assumed star formation history. The code also traces 
the energy and mass ejected from the stars through stellar winds 
and supernova explosions, as well as the UV radiation they emit. 
These quantities are the primary tools necessary to study the 
dynamical effects of the interaction between massive stars and 
the surrounding ISM. We interpolate between the discrete inputs 
from the models, mainly using linear interpolation, except for 
temperatures, masses, luminosities, surface gravities and stellar 
ages, which we interpolate logarithmically, unless otherwise 
noted. 

The code is divided into two separate parts: in the first one 
stellar tracks are used to calculate time profiles for individual 
stars, and we interpolate between these to create a fine grid 
of stellar tracks as a function of progenitor mass. For a given 
time we tabulate the stellar properties as a function of initial 
stellar mass, thereby creating isochrones. The second part of 
the code uses the isochrones as input, and for each isochrone 
calculates integrated quantities for a given population of stars, 
weighting the isochrone values by the IMF. This part of the pro- 
gram calculates analytically the mean value of the distribution 
of properties we are interested in, as well as some high order 



moments that allows an analytical estimation of the probabil- 
ity distribution of populatio ns with a finite (and small) size (see 
ICervino & Luridianal 120061, for details). It also estimates these 
distributions with Monte Carlo simulations for populations of 
any size, where intial stellar masses are chosen randomly, with 
relative probabilities according to the assumed IMF. These out- 
puts are complementary in the sense that the average values are 
suitable for very large populations of stars, such as entire galax- 
ies, whereas the Monte Carlo simulations are necessary for very 
small star-forming regions. The analytical probability distribu- 
tions are useful in intermediate cases ranging from about one 
hundred to thousands of massive stars. As we aim to model star 
forming regions in the solar neighbourhood, we restricted our 
study to stellar models of solar metallicity. 

The details of the physical ingredients of the code are de- 
scribed in the following. For a number of processes, switches are 
included to select between alternative models. We identify these 
through convenient parameter names that are used throughout 
the paper when discussing the outputs of the code. The paramer 
names are summarized in Table Q] 

2.1. Stellar evolution 

For the stellar evolution of massive stars several alternative sets 
of stellar tracks are implemented. The default (geneva05) con- 
sists of the solar metallici t y, rotating stell a r mod els described in 
IMevnet & Maederl (I2005I) : IPalacios et all d2005l) . These models 
all have ZAMS rotation velocities of 300 km s , producing 
time averaged equatorial velocities on the main sequence 
between 200 and 250 km s _1 . There are two ve rsions of these 
mode ls, one including the calculation of 26 Al (Palacio s et al.l 
I2005I) for stars with initial masses betwe en 25 and 120 M Q , 
and one without (Mev net&Maederl [2005) for stars with initial 
masses b etween 9 and 1 20 M Q . We therefore use the stellar 
tracks of Palacios et al. (2005h above 25M , combined with 
the models of IMevnet & Mae der (2005) below this limit as our 
default. However, to enable the calculation of 26 Al, parts of 
the numerical calculation was changed between the two sets 
of models (G. Meynet, private communication). We therefore 
test our results on the stellar energy ejection against those 
obtained using stellar tracks of Mevnet & Maeder (2005) only 
(geneva®5alt), see section [331 

The main effects of rotation on the yiel ds of 26 Al, and the 
energy and mass ejection a re the following (IMevnet & Maederl 
2005; IPalacios et all 120051) : Rotational mixing allows surface 
enrichment in 26 Al at an earlier evolutionary stage than obtained 
in models without rotational mixing. When rotation is not 
accounted for, surface 26 Al enrichment occurs only when the 
deep layers of the stars where 26 Al is synthesized are uncovered 
by the stellar winds. In rotating models rotational diffusion 
enables surface enrichments well before the regions processed 
by hydrogen burning are uncovered by the stellar winds. This 
effect increases the quantity of ejected 26 Al with respect to 
non-rotating models. For the same reasons rotating stars present 
surface abundance characteristics of the Wolf-Rayet stages 
before the stellar winds have uncovered the core. This causes 
the stars to enter the WR-phase earlier and to increase the WR 
lifetimes. For example a 40 (85) M e non-rotating star has a WR 
lifetime of merely ~0.1 Myr, whereas the rotating counterpart 
has a WR lifetime of almost 0.4 (1.4) Myr. Also rotational 
mixing enables stars with masses as low as 22M to become 
WR-stars. These effects allow rotating models to eject more 
26 Al into the ISM than non-rotating ones. This also increases 
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Parameter name 


Parameter description 


References 


genevaQ5 


stellar evolution with rotation 


Mevnet & Maeder (2005); Palacios et al. (2005) 


genevaQ5alt 


stellar evolution with rotation 


Mevnet & Maeder (2005) 


geneva97 


stellar evolution 


Maeder & Mevnet (1994); Mevnet et al. (1997) 


LC06 


stellar evolution 


Limongi & Chieffi (2006) 


windQ8 


wind velocities 


Lamers et al. (1995); Niedzielski & Skorzvnski (2002) 


windQQ 


wind velocities 


Howarth & Prinia (1989): Prinia et al. (1990) 


yieldsLCQ6 


supernova nucleosynthesis yields 


Limongi & Chieffi (2006) 


yieldsWW95 


supernova nucleosynthesis yields 


Wooslev et al. (1995); Cervino et al. (2000a) 


atmosMS 


stellar atmospheres 


Kurucz (1992): Martins et al. (2005): Smith et al. (2005) 


atmos Smith 


stellar atmospheres 


Kurucz (1992); Smith et al. (2005) 


atmosOLD 


stellar atmospheres 


Kurucz (1992); Schmutz et al. (1992); Schaerer & de Koter (1997) 



the mass loss and therefore the amount of kinetic energy ejected 
into the ISM. Furthermore stellar lifetimes are increased by 
rotation, with a 15-25% increase in the hydrogen burning 
lifetimes. This delays the onset of the supernova explosions, 
and lowers the supernova rate for a population of rotating stars, 
compared to a population of non-rotating stars. 

Rotational mixing can in many respects improve the cor- 
respondence between the outputs of stellar models and the 
observed features of massive stars. For instance, rotational 
mixing may be the cause for s urface nitrogen enrichments in OB 
main sequence stars (see e.g. iHunter et alj 120071; iMaeder et al.l 
2008). With the relatively low mass loss rates presently favoured 
(see below), non-rotating models underpredict the observed 
ratio of WR to O-type stars. Rotational mixing favors WR 
formation, and rotating models are thus in better agreement 
with observations. Provided that most of the nitrogen enriched 
stars and WR stars are produced by single star evolution, 
rotation appears as a key physical ingredient of the models. 
However, we note that the comparison of models with obser- 
vations and the discussion about whether the rotating models 
reproduce the observed properties of massive stars better than 
non-rotating model s is currently a subject of debate (e.g. 
Hamann et aTll2006t IMevnet & Maede rl 120051: lEld ridge & Vinkl 
20061; I Vazquez & Leithererll2007l) . 

Two alternative sets of non-rotating models are imple- 
mented. O ne set (geneya97) c onsists of the solar metallicity 
models of IMevnet et alj dl997b combined with the models of 
IMaeder & Mevnet! d 19941) for stars b elow 25 Af B . The other se t 
(LC06) consists of the stellar tracks o flLimongi & Chie ffi (2006). 
The main differences between the sets of stellar evolution mod- 
els relevant to our study are the inclusion of rotation in the 
geneva05 stellar models, and the updated wind mass-loss es- 
timates included in the geneva05 and LC06 calculations. These 
new estimates reduce the wind mass-loss rates of the geneva®5 
and the LC06 models significantly compared to the geneva97 
models: in the Wolf-Rayet (WR) phases the inclusion of clump- 
ing in the interpretation of the observational data leads to re- 
duced mas s-loss rates by a fact or of 2 - 3 compared to the earlier 
models dNugis & Lamer s 2000). Als o the theore tical mass-loss 
rates for O and B stars of lvink et all d2000l 1200 lb are lower by 
similar factors in co mparison to earlier empir ical values of e.g . 
Ide Jager et ail (Il988l) . See the recent review bv lPuls et all d2008b 
for an extensive description of mass-loss rates. 
We follow each stellar track through the provided evolutionary 
points and calculate the quantities needed for the population syn- 
thesis. 



2.1 .1 . Kinetic energy 

The rate of kinetic energy emitted in the winds of the stars is 
calculated from E^-l/lMv^, where M is the mass-loss rate 
from the stellar atmosphere and is the velocity of the wind 
at infinity (the terminal velocity). This represents the energy 
available once the wind has escaped the gravitational potential 
of the star. The interaction of the winds with the ISM can 
convert a large fraction of the energy into gas turbulence and 
radiation. The mass loss rate is given by the stellar tracks. To 
calculate the wind velocity we firs t coarsely classify th e stars 
according to the following criteria (Leithe rer et al.lfl999h : Stars 
with a mass loss rate above M = 10~ 3 ' 5 M o yr _1 and effective 
temperatures in the range 3.75 > logr e jf < 4.4 are classified 
as luminous blue variables (LBVs). Stars with an effective 
temperature above log T e ff = 4.0 and a fractional abundance of 
hydrogen at the surface below 0.4, are considered WR-stars. 
They are furthermore divided into subc l asses a ccording to the 
surfac e abundances ( Smith T& Maederl |199 it iLeitherer et al.l 
1 1999b . see Table |2] 

The terminal velocity of the wind depends on the type of 
star. We utilize two different prescriptions for the calculation: 
in the default mode (wind®8) stars outside the categories de- 
fined above are divided into hot stars with v M = 2.6v e sc and 
cool stars with v^, = 1.3v esc at a temperature of log T e ff = 4.32 
(Lamers et al. 1995), where v esc is the escape velocity at the sur- 
face of the star. Luminous blue variables are assumed to have 
winds with low velocities (20 km s" 1 ) and the velocities of the 
WR- winds are estimated from Niedzielski & Skorzvnskil (f2002) 
and given in Table [2] 

The alternative mode windSS corresp onds to the assumptions in 
Starburst99 dLeitherer et al.ll 1992b . Here stars outside the cat- 
egories have w ind velocities given by V oo = 0.85 x (0.58 + 2.04 x 
log x Vesc dHowarth & Prinia1ll989h . luminous blue variables 

again have = 20 km s" 1 , and the v elocities of the WR-winds 
are estimated from Pri nia et ail (Qj)90), and given in Table [2] 

2.1.2. Radioactive isotopes 

The wind ejection of the radioactive isotope 26 Al is followed 
along the stellar tracks. The ejection rate is found from multi- 
plying the surface abundance fraction of 26 Al with the mass loss 
rate, and the amount of 26 Al present in the interstellar medium 
around the star is calculated by taking into account the decay 
timescale of ~ 1 Myr. The other radioactive isotope considered 
in this study, 60 Fe, is not ejected in the stellar winds, but only 
in the supernova explosions. For this isotope we consider two 
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Table 2. Classification criteria for WR stars, and their wind velocities. *H S ,C S ,N S and He s are the fractional surface abundances 
(mass fraction) of hydrogen, carbon, nitrogen and helium, respectively, and the wind velocities are in km s _1 . 



WR type 


Surface abundances 




Wind velocity (wind08) 


Wind velocity (windSO) 


WNL 


0.4 > H s > 0.1 




1250 


1900 


WNE 


H s < O.l&C/W, < 10 




2000 


1650 


WC6-9 


H s < 0.1&C S /Ns > 10&(C S + O s )/He s 


<0.5 


1760 


1810 


WC4-5 


H s < 0.1&C S JN S > 10&(C S + O s )/He s 


< 1.0 


2650 


2820 


WO 


H s < 0.1&CJN S > 10&(C S + O s )IHe s 


> 1.0 


3000 


3000 



decay timescales, with ~2 Myr as the default and ~3.8 Myr sug- 
gested by recent measurements. At the last point of the stellar 
tracks the radioactive yields from the supernova explosions are 
added. We have two prescriptions for calculating the yields: in 
th e default mode (yieldsLC 2686) we use the yields calculated 
by lLimongi & Chieffil (12006b . The yields are found based on the 
initial mass of the stars, interpolating linearly between their data 
points. In the alternative method (yields WW95) the yields be- 
low in itial masses of 25M are taken from [Wooslev & Weaver] 
(119951) . Above this mass we set the output of each star to 
5.0 • 10" 5 M o for both 2 6 Al and 60 Fe. This is a pproximately the 
average yields found in ICervino et alj d2000al) from linking the 
core mas s of the geneva97 s tellar tracks with the supernova 
yields o f IWooslev et al.l ( 1 1995b . We note that in ICervino et alj 
(2000a), the results depend strongly on the method used (see 
their Figs. 1 and 2). Our crude approximation therefore reflects 
unknown physics. The total yield is not affected dramatically by 
these uncertainties, as there are few massive stars, and for 26 Al, 
the emission from these stars is dominated by the wind contri- 
bution. 

2.2. Creation of the stellar isochrones 

We use the stellar tracks to create isochrones for a user-defined 
set of times after the onset of star formation. To facilitate the 
interpolation between the stellar tracks, for each of the tracks a 
subset of 5 1 equivalent evoluti onary points ar e extracted accord- 
ing to the definitions given in Maeder (1990). The reduction of 
the stellar tracks to 5 1 points eliminates some of the finer details 
of stellar evolution, but the uncertainties due to this simplifica- 
tion are small compared to the uncertainties of defining equiva- 
lent evolutionary points and interpolating between stellar models 
with different initial mass. For each isochrone an array of initial 
stellar masses is defined, and stellar tracks are produced by inter- 
polation between the equivalent points of the nearest input stellar 
tracks. We note that the ages at the equivalent points are interpo- 
lated using a cubic spline in logarithmic space, as lower degree 
interpolation methods create artif icial results (e.g. in the super- 
nova rate, see ICervino et alJfeOOll Fig. 1). The isochrone values 
are then found by interpolation along these stellar tracks. We 
note that the geneva05 tracks are only calculated to the end of 
central helium burning, while the geneva97 tracks are followed 
to the end of central carbon burning. We use the last available 
grid points and extrapolate it towards the end, thus neglecting 
the stages between the ends of the stellar tracks and the super- 
nova explosions. 

2.3. Calculation of population properties 

The stellar isochrones are the basis of our population synthesis. 
Each isochrone is used to calculate the population properties, 
both by integrating the isochrone analytically and using Monte 



Carlo techniques. Both methods sample the number of stars of 
a given mass according to the initial mass f unction (IMF), for 
which we use the Salpeter mass function (|Salpetenfl955h as 
the default, with 8M and 12OM as lower and upper integra- 
tion limits. As only massive stars are investigated here, we do 
not consider more complicated IMF models (e.g. IScaiol ll986: 
lKroupal2001h that also address the shape of the mass function at 
lower masses. We furthermore i nclude an implementation of the 
statistical methods described in ICervino & Lurid iana (200g) to 
account for statistical variability of stellar-mass weightings for 
smaller numbers of stars. This allows us to calculate analytically 
the statistical distribution of the output quantities, depending on 
the number of stars in the stellar population, similar to what can 
be done with Mont e Carlo simulations. We do not include th e 
effect of binaries dVanbeveren et alj|2007t lEldridge et afl l2008). 
The frequency of interacting binaries, and the effects of the inter- 
actions are poorly known, and the effects on the parameters stud- 
ied in this paper are difficult to assess. Such interactions might 
remove outer layers of primary stars, and rejuvenate secondary 
stars. Both process es can potentially enhance the 26 Al ejection 
dLangeretalJll998h . 



2.3.1. Supernova explosions 



In this study, we assume that the kineti c energy ejected in a 
supernova explosion is 10 51 erg s _1 (e.g. Woo slev et alJll99"5h . 
irrespective of the progenitor mass. The supernova rate is 
derived as an average over a period of 0. 1 Myr and from this the 
'kinetic luminosity' of the supernovae is calculated. The mass 
ejected in supernova explosions is found from the difference 
between the mass of the last point in the stellar evolution code, 
and the mass of the compact remnant. We assume that all 
neutron stars have a canonical mass of 1 .4M and that all black 
holes have a mass of 7.OM . The supernova explosions of stars 
with initial masses below 25M are assumed to produce neutron 
stars whereas stars above this mass limit are assumed to produce 
black holes. Our results are not sensitive to these values within 
reasonable ranges. For example, if the most massive stars end up 
as neutron stars due to mass loss, this will only change the mass 
of the ejected matter by ~ 5 - 10%, as the majority of the mass 
is ejected in the winds at earlier evolutionary stages. It should 
be noted that the geneva97 and geneva®5 stellar tracks do not 
cover the final stages of stellar evolution, and the mass ejected 
in the supernova explosions may therefore be overestimated. 
However, the mass loss from stellar winds is underpredicted by 
the same amount, and these effects should cancel out to some 
extent. 
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Time [Myr] 

Fig. 1. Comparison between Starburst99 (thin black lines) 
and our results (thick grey lines) for a co-eval population, top: 
time profiles of the kinetic power emitted through the stellar 
winds (solid) and supernova explosions (dotted). Bottom: The 
ionizing flux. 



2.3.2. Ionizing flux 

The ionizing flux of each star is calculated by linking stellar at- 
mosphere models to the state of the star. The surface gravity 
and the effective temperature are used to find the atmosphere 
model with th e closest res e mblan ce. In the WR-phases, we use 
the method of ISmith et alj d2005l) to find these values at an op- 
tical depth of r = 10. In this method, the effective temperature 
is simply given by T» = 0.6 ■ 7h y d + 0.4 • T2/3, where 7h y d is 
the hydrostatic effective temperature of the stellar models, and 
T2/3 is the corrected surface temperature given in the geneva97 
and geneva® 5 stellar models. This parameter is unavailable in 
the L&CQ5 stellar models, and the ionizing flux is therefore not 
calculated for these. 

The defau lt set of atmospher e models (atmosMS) applied for 
OB s tars dMartins et al.l 120051) and for WR stars dSmith et alj 
120051) both take non-LTE, wind and line-blanketing effects into 
account. For sta rs with param eters beyond the range of these 
models we use iKuruczl d 19921) atmospheres. For comparison 
wit h Starburst99 we include an option of using the models 
of dSmith et al.ll2005l) for both WR and OB stars (atmosSmith), 
and for comparison with earlier studies an alternat ive set of mod- 
els (a tmosOLD) consisting o f the WR models of Schmutz et alj 
d 19921) and the OB models of lSchaerer & de Koterl d 19971) . again 
supplemented by the models of lKuruczl d 1992b . 



3. Population synthesis results 

We combine the results obtained from integrating over the indi- 
vidual isochrones to create time profiles of the average ejection 
of 26 Al, 6l) Fe, kinetic energy, hot gas and the ionizing flux. In this 
section we discuss these individual quantities, and their depen- 
dence on the assumptions made. 



3. 1 . Comparison with previous studies 

Synthesis of stellar populations is an important tool for un- 
derstanding stellar clusters and galaxies. Since the first studies 



(lTinslevlll978t lBruzuallll983l) t he sophistication of the m odels 
have increased drastically (e.g. iBruzual & Chariot! I2003). The 
majority of these models are designed to predict the spectra of 
large stellar populations over a wide range of ages and metallic- 
ities, for which the exact parameters of the few massive stars 
are not of major importance. This is quite different from the 
goal of our study, which is the calculation of parameters of in- 
terest for the interactions between massive stars and their envi- 
ronment. Most studies are therefore not suited for comparison 
with the present study. Thus we use Starburst99 for compar- 
ing our results on the ionizing flux and kinetic luminosity. For 
this purpose we produce a model with parameters similar to pa- 
ramet ers of their stu dies. We choose a model with a Salpeter 
IMF (Sal peter! [T955I) . between 8 - 12OM with the parameters 
geneva9/Q, wind®®, atmosSmith and yieldsWW95. The com- 
parison for a cluster with a total mass of 1O 6 M in the given mass 
range are shown in Fig. Q] There is obviously good agreement 
between both models. The few discrepancies (e.g. the somewhat 
larger wind power between 2 and 3 Myr in our model) can all 
be traced back to different interpolation schemes, and therefore 
represent unknown physics due to the limited stellar grid. The 
o nly previously publishe d population synthesis of 26 Al and 60 Fe 
is ICervino et al.l (l2000al) . Again we find good agreement with 
their results. Only the early supernova ejection of 26 Al is higher 
in our model. This is due to the fact that their supernova yields 
decrease for the highest mass stars, whereas for our parameter 
yieldsWW95 they are constant above 25M . 

3.2. Radioactive isotopes 

Compare d to the previous po pulation synthesis study of 26 Al 
and 60 Fe dCervifio et alj |2000a), our study includes new stellar 
models, specifically geneva® 5 and LCQ6. Both take the reduced 
mass loss estimates into account, and the geneva® 5 model fur- 
thermore includes the effects of stellar rotation. The LCQ6 model 
consistently includes the supernova explosions of the very mas- 
sive stars. While these improvements have been discussed for 
individua l stars and for the overall ejection from a p opulation 
of stars dPalacios et alj|2005t iLimongi & Chieffill2006l) . the ef- 
fects on the time profile have not been evaluated in detail. The 
time profiles of 26 Al and 60 Fe in the interstellar medium for a co- 
eval population of stars are shown in Fig. [2] The 26 Al has been 
divided into the wind and supernova contributions. Only the su- 
pernova contribution to 60 Fe is shown, as the wind contribution 
is negligible. 

The wind contribution of 26 Al is relatively similar for the 
geneva®5 and the geneva97 models. This is quite surprising, 
given the relatively large differences between these models. The 
reason is that the lower wind mass-loss in the geneva® 5 models 
is compensated by two effects. One is the rotational mixing of 
elements, which leads to larger surface abundances. The other 
is the longer lifetime of the rotating models, which causes the 
amount of material lost to be similar to the amount lost in the 
geneva97 model, despite the lowerrate of mass loss. This is also 
the reason that the 26 Al time profile shows a more extended peak 
for the rotating models. The LCQ6 models, on the other hand, 
predict much less 26 Al ejection by the pre-supernova wind, for 
two main reasons: with the reduced mass loss, the rate of ejec- 
tion is low, and the 26 Al is not effectively mixed to the surface. 
For the supernova ejections the differences between the model 



1 While t he genevaQ5alt stellar models are included in 
Starburst99 (Vazquez & Leitherer 2007), this part of the code is not 
publicly available yet. 
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Fig. 2. Average time profiles of 26 Al and 60 Fe for a coeval popu- 
lation of stars, for the three different sets of stellar models avail- 
able. The lines show the amount of the elements present in the 
ISM per star in the 8 - 12OM range. The 26 Al yields are divided 
between the wind and the supernova contributions, whereas the 
wind contribution of 6() Fe is negligible. The black lines indicate 
the results using the yieldsLC2006 supernova yields, whereas 
the grey lines indicates the results using the yieldsWW95 yields. 



predictions of Al and Fe are smaller, with the main effect 
coming from the longer lifetime of the rotating models causing 
a slightly less peaked profile. On the other hand, drastic differ- 
ences between the predictions for the first 10 Myr are obvious 
when comparing yieldsWW95 and yieldsLC2006. This proba- 
bly mainly reflects the uncertainties in the final structure of the 
very massive stars, at the time o f the supernova explosion s. The 
low mass-loss in the models of dlimongi & Chieffi| [2006) gives 
very massive cores, that produce large amounts o f the elements. 
The h igher mass-loss rates assumed in the study of lCervino et al.l 
(2000a|) gives lighter cores, and this results in the yieldsWW95 
being lower than the yieldsLC2Q06. We note that with the very 
large theoretical uncertainties, observations of 26 Al and 60 Fe thus 
have the potential to place interesting constraints on the final 
evolutionary state of very heavy stars, although in an indirect 
way, since only the integrated effect of different sources can be 
observed. 

In Fig. [3] the average time profiles of 26 Al (including both the 
wind and the supernova contributions) and 60 Fe are shown for 
different models, together with the statistical variance. Shown 
are the lcr and 2<x statistical deviations (the intervals containing 
68% and 95% of the Monte Carlo simulations) for a population 
of 100 stars in the 8 - 120M o rangfl corresponding to a typical 



2 1 star in this range corresponds to 381, 140 an d 188 sta rs in the 
0.1 - 120M o range for the Salpe ter dSalpeteilll955h . Kroupa dKroupal 
1200 lh and the Scalo (Scalo 1986) mass functions, respectively, and to 
13%, 18% and 5% of the stellar mass. 



nearby star forming region (for example the number of massive 
stars formed within the last 15 Myr in the Orion O B 1 association 
is estimated to be close to lOO lBrown et al.|[l994l) . We note that 
the 8-120M Q range includes more massive stars than are ob- 
served in many nearby regions. However, when a probabilistic 
description (such as our Monte Carlo simulations) is assumed, 
the limit should be the most massive star theoretically possible 
in the cluster. Observed clusters correspond to random realiza- 
tions of the IMF and the most massive stars in these can therefore 
have much lower masses than the upper limit of 12OM (see also 
section [4] where the contributions of various initial mass ranges 
to the observables are shown). 

From Fig. [3] it is obvious that for relatively small populations, 
it is essential that these statistical effects are taken into account, 
when interpreting observations. Also very interesting is the ra- 
tio between the observable strengths of the 6l) Fe and 26 Al y-ray 
lines. This is shown in Fig.|4]for the same models as in Fig. [3] 
For much of the time, this ratio places stronger constraints on the 
stellar models than the individual observations of 26 Al and 60 Fe. 
This is due to the fact that the emissions of these two elements 
are correlated. The strong increase in the 60 Fe/ 26 Al ratio seen 
around the lifetime of an 8M star (~ 35 - 50 Myr, depending 
on the stellar model) is simply an effect of the longer lifetime of 
60 Fe, when the elements are not being replenished (non-steady 
state situation). 

Recent results (G. Rugel, in preparation) indicate that the life- 
time of 6l) Fe is significantly (~3.8 Myr) longer than the com- 
monly used lifetime of ~ 2 Myr. It is unclear if a different life- 
time would have any significant impact on the nucleosynthesis 
of Fe, and a study of this is beyond the scope of this paper. 
In Fig. [5] we show the effect on the time profile of 60 Fe for 
our default model, assuming that the amount of 60 Fe released 
in the supernova explosions is unchanged. Due to the longer de- 
cay timescale, the build-up of the isotope is larger, and in the 
period after the main peak (5-15 Myr after the star forma- 
tion), the amount of 60 Fe present in the ISM is approximately 
twice as large as for the shorter 60 Fe lifetime. Note that the ef- 
fect on the observed flux is different. While there is more 60 Fe 
present in the ISM, the y-ray emission per unit mass is decreased 
due to the longer decay timescale. This is illustrated by the grey 
dashed line in Fig.|5j the integrated amount of y-ray emission is 
unchanged (since the amount of 60 Fe released from the stars is 
unchanged), but the distribution is slightly shifted to later times. 
The only significant effect is the lower peak at ~ 5 Myr. It should 
be noted that in a constant star-formation scenario (steady-state), 
the measured ratio of 6l) Fe/ 26 Al emission is independent of the 
decay times of the two isotopes (although the ratio of 60 Fe/ 26 Al 
does depend on the decay times). 

3.3. Energy and mass 

The energy and mass ejection rates are the most important 
parameters for the evolution of the interstellar medium in the 
vicinity of young stars. Winds and SN ejecta compress the 
surrounding medium, creating wind-driven bubbles, supernova 
remnants, and even super-bubbles in the case of multiplets of 
massive stars. These stellar systems dissociate the surrounding 
molecular clouds and may induce star-formation within these. In 
Fig. [6] the "kinetic luminosities" of the winds and the supernova 
explosions are shown for the three sets of stellar tracks. In 
the first 3-4 Myr, no stars explode, and the kinetic luminosity 
originates exclusively from stellar winds. The mechanicial 
power of these sustained winds is roughly equal to the average 
power of punctuated supernova explosions, once these set in. 
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Fig. 3. Time profiles of 26 Al (top) and 60 Fe (bottom) for a coeval 
population of stars. The solid lines indicate the average profiles 
for the geneva05 stellar models with the yieldsLC2O06 su- 
pernova yields. The the dark and light grey regions show the lcr 
and the 2<x deviations for a population of 100 stars between 8 and 
120 M Q , based on Monte Carlo simulations. The dashed and dot- 
ted lines show our main alternative models. Dashed: geneva97 
stellar tracks with yieldsWW95 supernova yields. Dotted: LC06 
stellar tracks with yieldsLC2(S(56 supernova yields. 
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Fig. 4. Time profiles of 60 Fe/ 26 Al emission ratio for a coeval pop- 
ulation of stars. Legends same as Fig. [3] 
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Fig. 5. Time profiles of 6l) Fe for a coeval population of stars for 
the two different decay timescales. The solid line is the same as 
in Fig. [3] with an average lifetime of ~ 2 Myr, whereas the black 
dashed line is for an average lifetime of ~3.8 Myr. The grey 
dashed line indicates the observational difference (the amount of 
60 Fe that would be inferred from observations if a lifetime of 2 
Myr is assumed, in the case the lifetime actually is 3.8 Myr. 



ferent types of massive stars still have very large uncertainties, 
and this translates into large uncertainties in the implicated 
wind power. As the most massive stars explode after ~ 5 Myr, 
the wind power decreases rapidly, thereafter the supernova rate 
decreases steadily until stars with masses ~ 8M explode, and 
the energy deposition in the ISM becomes negligible. 

In Fig. [7] the average total kinetic luminosity and the mass 
ejection rate are compared through Monte Carlo simulations of 
a small population of stars. For the Monte Carlo simulations the 
power was averaged over 1 Myr, and still the variations are very 
large, especially in the later stages, where only few supernovae 
are likely to explode within a time interval. When studying small 
regions it is important to note these statistical effects, especially 
for processes with shorter timescales, where the variability of a 
local supernova rate increases drastically. For example the lcr 
range of the supernova rate around 10 Myr covers more than a 
factor of 100 when averaged over a timescale of 0. 1 Myr, similar 
to the replenishment times cale of the X-ray e mitting gas in the 
region of Orion studied bv lGiidel et al.l d2008l) . W e note that the 
variation i n the time integrated energy (studied bv lCervino et al.l 
1200 lL [2002b . important for studying the X-ray bubble expansion 
at larger scales, is smaller, since it only depends on the number 
of supernova explosions and not on the specific explosion times. 
The mass ejection rate shown in the lower panel of Fig.|7]is more 
strongly dominated by the ejection through stellar winds, and as 
this is a continuous process, the variations are relatively smaller. 



It is important to note that integrated over the first 10 Myr, the 
total energy emitted in winds is actually larger than the energy 
associated with the supernova explosions. This is important 
for shaping the medium around a star-forming region. There 
are large differences between the various stellar models. In 
particular it is worth noting that while decreased mass-loss rates 
in rotating stellar models obviously decrease the wind power, 
this is compensated for by their longer lifetimes. Furthermore 
the longer stellar lifetimes lower the supernova rate (except at 
late times >35 Myr, where the supernova rate for non-rotating 
stars becomes zero, whereas the rotating stars still explode until 
~45 Myr), and therefore the early dominance of wind energy is 
even stronger. 

A comparison between the two wind prescriptions windQQ and 
wind08 shows significant differences. Wind velocities of dif- 



3.4. UV radiation 

The bulk of the ionizing flux is provided by hot stars that disap- 
pear after only a few Myr. This is clearly seen in Fig. [8] where 
the time profile of the ionizing flux is shown. The flux is high 
and roughly constant during the first 2 Myr, after which it de- 
clines. After less than 10 Myr it has declined by a factor of more 
than 100, and the rate of decline steepens further. There are still 
large uncertainties in the modelling of non-LTE expanding at- 
mospheres. Comparing the results with a now superseded set of 
stellar atmospheres shows only minor differences, but the com- 
parison with the non-rotating L&C stellar models shows dramatic 
differences. So while the correct atmosphere modelling is impor- 
tant for the study of single stars and for understanding the UV 
spectrum, the stellar evolution modelling is much more impor- 
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Fig. 6. Time profiles of the wind (top) and supernova power 
per star from a coeval population of stars in the 8 - 12OM 
mass range, using the various stellar models with the wind08 
wind velocities. The solid dark grey line is associated with the 
geneva05alt evolutionary tracks, and the solid light grey line 
corresponds to the geneva05 stellar models with windOO. 



tant for the understanding of the interactions between groups of 
massive stars and their surroundings. It is important to note that 
for small stellar associations, the statis tical deviations from the 
mean are very large (see e.g. Fig. 4 in ICervifio et al.ll2000bl for 
the deviation for associations of different masses). 

3.5. Statistical properties of the stellar populations 

Above we described the average outputs of the 26 Al, 60 Fe, en- 
ergy, mass and UV radiation from young stellar populations, to- 
gether with examples of their variability, based on Monte Carlo 
simulations. For a better understanding of the statistical proper- 
ties of the s tellar distributions, we also m ade use of the statistical 
methods of ICervifio & L uridiana (2006). At a given time, these 
methods allow the analytical approximation of the probability 
distribution of a quantity. We calculate the raw moments of the 
distribution 



Time [Myr] 



Fig. 7. Time profiles energy (top) and mass (bottom) ejection 
rates, with lcr (dark grey) and 2<x (light grey) Monte Carlo devi- 
ations for a population of 100 stars between 8 and 120 M Q . The 
solid line indicates the full outputs, with the dashed and dotted 
lines indicating the contributions from stellar winds and super- 
nova explosions, respectively. The Monte Carlo outputs were av- 
eraged over 1 Myr. 
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where n is the integer order of the moment, I is the physical 
quantity for which we investigate the distribution, and <p(l) is 
the probability for a specific value of /. The parameter a is the 
average value of I: 



Fig. 8. Time profiles of the ionizing flux per star, with lcr (dark 
grey) and 2cr (light grey) Monte Carlo deviations for a popula- 
tion of 100 stars between 8 and 120 M , using the geneva05 
stellar models with the atmosMS atmosphere models. The thin 
solid line was produced by using the atmosOLD atmosphere 
models, whereas the dashed line uses the LCQ6 stellar models 
and the atmosMS atmospheres. 
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From these we find the skewness 



n = 



v 2 



3/2 



(2) 



(3) 



and the kurtosis 
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of the distribution. For an ensemble of stars, the distribu- 
tion expressed in normal form (i.e., transformed to a distribution 
with zero mean and unit variance) has a skewness Ti - j\j "V^V 
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and kurtosis F2 = ji/N. With these we can estimate the distri- 
butio n using a Gaussian multiplied by Edgeworth's series (se e 
e.g. iBlinnikov & Moessnerl 1 1998c ICervino & Luridia na 2006), 
which, truncated to 2 terms, is: 

(p Uot (x) = ——e x 
V27r 

1 1 + \ Tj (x 3 -3x) + ^- F 2 (x 4 - 6x 2 + 3) + 
\ 6 24 

J_r 2 (x 6 -15x 4 +45x 2 -15)). (5) 

For very small populations the discrepancy between the 
Edge worth approximation and th e real distribution becomes 
large. ICervifio & Luridiana (2006) therefore derived analytical 
estimates of the necessary population size to allow the use 
of their approximation, and also when it is possible to use a 
Gaussian approximation. In Fig. [9] an example of the probabil- 
ity distribution for 26 Al is shown, and the Edgeworth approxi- 
mation is compared to Monte Carlo simulations. It can be seen 
that while for very low numbers of stars, the distribution can 
only be investigated using Monte Carlo simulations, whereas 
for intermediate cases Edgeworth's approximation is useful. Fig. 
[TOl shows the skewness and the kurtosis of the distribution. In 
Fig. QT|the minimum number of stars needed for an adequate 
approximation (within 10%) in the +3<x interval is shown for 
Egdeworths approximation and for the Gaussian approximation. 
Evidently star-forming regions need to have a very large num- 
ber of massive stars (in the best case 1000 stars above 8 M ) 
for Gaussian statistics to be applicable. But even outside the 
Gaussian regime, the evolution of y\ and 72 values allow us 
to obtain some inferences about the results of the Monte Carlo 
simulations with low numbers of stars without performing them. 
First of all, j\ is a measure of the asymmetry of the distribution. 
A positive j\ means a "L-shape" distribution. A large j2 means 
that the distribution is both more peaked than Gaussian distribu- 
tions near the maximum, and more flat than Gaussian in the tails. 
Large j\ and j2 values together reflects then (a) the most prob- 
able value is peaked in a value different from the mean and, (b) 
there is a non-negligible high luminosity tail in the distribution. 

Appling these concepts to our study, the large values in y\ 
and 72 at young ages (before the onset of the SN phase in 26 Al 
production) shows the large impact of particular high massive 
stars in the the integrated production of 26 Al. It implies that, for 
low mass clusters, the current 26 Al produced by an individual 
cluster is strongly dependent on the most massive star the clus- 
ter has formed, so (provided the cluster age is known from an in- 
dependent method) it would be possible to estimate the value of 
this most massive star. Of course, since we are only able to obtain 
statistical inferences in the low cluster mass regime, we need to 
take into account the maximum possible M up value, which nat- 
urally include the case where such more massive stars are not 
present in the cluster due to the particular realization of the IMF 
performed by Nature in a given cluster. These effects are clearly 
visible even for cluster with 10 massive stars in the top pannel 
of Fig. 9. In the pure SN dominated phase (ages larger than 20 
Myr), j\ and 72 are almost constant, without significant structure 
and are both slowly increasing with age (cf. Fig fTTb . A simple 
consequence of the slow decay of the SN rate at these ages is, 
(a) an increase of sampling effects due to the declining number 
of possible SNe within a given time window, and (b) a decrease 
in the sensitivity to which star has exploded as a SN in the 26 Al 
ejection, as there is not much variation in the 26 Al ejection from 
SNe at these late times. The intermediate age range between 7 
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Fig. 9. The statistical distribution of Al in the interstellar 
medium around a cluster of 10, 50 and 200 massive stars (A'*) 
in the 8 - 12OM range in the top, middle and bottom panels, re- 
spectively. The distribution is given for an age of 6 Myr. The his- 
tograms show the distribution of 2000 Monte Carlo simulations, 
while the solid lines are based on the Edgeworth approximation. 
The vertical dashed line is the average value. geneva®5 stellar 
models and LC2006 supernova yields were assumed. 



an 20 Myr shows how the stars that explode in this age range 
have different amounts of 26 Al ejections. 

Finally, in Fig. [12] the Pearson correlation coefficient be- 
tween the y-ray emission from nucleosynthesis products and the 
hydrogen ionizing radiation is shown. The use of these correla- 
tion coefficients allows us to establish the probability of obtain- 
ing the value of one of the observables, given that the value of 
the other observable is known. A correlation coefficient equal 
to zero implies that the two associated quantities do not contain 
information about each other. There is no strong correlation be- 
tween the y-rays from 6() Fe and UV radiation, and therefore the 
overall UV radiation does not know anything about the 60 Fe pro- 
duction. The strong correlation between the 26 Al 7-ray emission 
(from the 26 Al ejected in the stellar winds) and the UV radiation 
at young ages refects that both components are produced by re- 
lated sources: Note that there is a time delay between the UV 
and the 26 Al emission, so the sources that produce both quan- 
tities are not the same. The 26 A1-UV correlation decreases with 
time, since the stars that produce the UV emission are unrelated 
to the stars that explode as SN and produce most of the 26 Al 
at ages later than 6 Myr. Finally, the y-rays from 60 Fe and 26 Al 
are highly correlated but the correlation does not reach a value 
of 1, reflecting the fact that the overall production comes from 
similar, but not the same, stellar sub-populations. 
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Fig. 10. The skewness, j\ and the kurtosis, 72 of the 26 Al distri- 
bution for a population of stars in the 8 - 12OM mass range. 
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Fig. 11. Lower limits on the number of stars in the 8 - 12OM 
range needed for the statistical distribution of 26 Al to be approx- 
imated by a Gaussian (upper line) and by Edgeworth's approx- 
imation (lower line). geneva®5 stellar models and LC2006 su- 
pernova yields were assumed. 
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Fig. 12. The Pearson correlation coefficient between the 26 Al and 
60 Fe y-ray emission (solid line), the 26 Al y-ray and the UV emis- 
sion (dashed line) and the 60 Fe y-ray and the UV emission (dot- 
ted line), as a function of time for our default model. 



studies, as they are a unique way of tracing the output from the 
massive stars directly, globally (Milky Way wide) as well as the 
spatial distribution of the ejecta around young stellar clusters. 
Many aspects of stellar evolution of massive stars are quite 
uncertain, and often the effects of these on the properties of a 
stellar cluster are not directly apparent. Population synthesis 
studies allow to test the impact of changes in specific ingredients 
and processes of stellar evolution models. 
Our study of the ejection of the radioactive elements 26 Al 
and 6(, Fe shows considerable theoretical uncertainties on the 
ejection rates and time profiles. Especially the amount of 60 Fe 
ejected from the most massive stars is sensitive to the structure 
of the stars in the final evolutionary stages. The emission of 
26 Al varies by less than a factor ~2, and the ratio of these 
two elements is therefore an important diagnostic for the late 
stages of stellar evolution of very massive stars. An important 
ingredient is mass loss from the massive stars. Despite recently 
improved understanding of this process, the rate of mass loss is 
still subject of debate, with proposed downward changes of up 
to an order of magnitude. One effect of lower mass loss rates 
is an increase of the stellar core sizes at late burning stages. 
The 6l) Fe production is correspondingly enhanced, and the wind 
ejection of 26 Al is decreased. For non-rotat ing stellar models, 
with the currently preferred mass-loss rates dLimongi & Chieml 
2006) the Galacti c 60 Fe/ 26 Al ratio i s overpre dicted with respect 
to me asurements (Wang et al. 2007), leading Limongi & ChiefR 
(2006) to suggest changes in the IMF or a relatively low upper 
integration limit. For lower mass-loss rates this problem would 
be enhanced and this is a strong argument against very low 
mass-loss rates for non-rotating stellar tracks (although it is 
important to note that the uncertainties in the nucleosynthesis 
of 6l) Fe and 26 Al are too large for the argument to be definitive 
at this time). Stellar rotation can strongly enhance the mass lost 
from stars, mainly from increasing the duration of the WR- 
phases. However, rotation also tends to increase the size of the 
stellar core, and without calculations of the later evolutionary 
stages of such stars, it is difficult to predict the explosively 
ejected nucleosynthetic yields. The observed average Galactic 
60 Fe/ 26 Al ratio provides interesting constraints on the models, 
but observations of specific regions are also important, as 
the time profile and the statistical variability can improve our 
understanding of massive star evolution considerably. While 
the study of the ratio in individual star-forming regions will not 
be possible for many years, it is expected that it will at least 
be possible to divide the Milky Way into a few sections for 
this purpose at the end of the INTEGRAL mission. The main 
problem is the detection limit of the 60 Fe y-rays. On the other 
hand, the prospects for studying the time profile of 26 Al are 
good, and in subsequent papers we will apply our population 
synthesis tool to the Galactic 26 Al distribution, and to individual 
star forming regions, such as Orion. Fig. [13] shows that at 
different intervals, the output of energy and matter from the 
stars is dominated by different ranges of initial stellar masses. 
Therefore the analysis of different regions with non-steady star 
formation can be used to explore the physics of stars in quite 
narrow mass ranges. 



4. Discussion 

The outflow (energy and matter) from young stars into the ISM 
determines the interplay between star formation and galaxy 
evolution. Our population synthesis code predicts important 
stellar outputs necessary to study these processes in more detail. 
The emission of radioactive isotopes is important for such 



Al and Fe, together with many additional isotopes, are 
emitted into the ISM with very high velocities. These ejecta will 
move away from the star clusters, which has been observed in 
the Orion region where a part of the 26 Al signal is off-set from 
the star clusters, and coinc ides with the Eridanus bubble formed 
by the outflowing hot gas (Diehl 2002). The spread of 26 Al can 
therefore trace the mixing of young stellar ejecta into the ISM, 
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a process that is relatively poorly understood but important for 
proper modelling of the evolution of galaxies. Our code predicts 
the mass, energy and ionizing radiation output of the stellar clus- 
ters, which are the most important parameters for the interaction 
with the ISM. Our results show that stellar winds are important 
for the turbulent state of the ISM. In the first ~6 Myr the power 
from the winds clearly dominates, being as strong as the subse- 
quent power injected by supernovae, even if the total energy of 
the supernova explosions is higher than the total wind energy, 
as they explode over a much longer timescale. For the winds 
the most massive stars are very important, even if there are rela- 
tively few of them, as they lose most of their mass through winds 
with very high velocities in the WR-phases. Estimating the wind 
power from a typical O-type star and multiplying by the num- 
ber of stars in a cluster will underestimate the wind power by 
more than an order of magnitude. When comparing the mechan- 
ically derived power from different stellar models, we find that, 
while the mass-loss rates are smaller by a factor ~ 2 - 3 in the 
genevaQ5 models compared to the geneva97 models, the wind 
power is actually more dominates with the geneva05 models. 
On the other hand, for the LC06 models that also have the lower 
mass-loss rates, the role of the winds is somewhat diminished. 
The reason for the enhanced wind importance when using the 
rotating geneva05 models is caused by several factors, the most 
important ones being the greatly increased time spent in WR- 
phases where the mass-loss rate and wind velocity is very high, 
and the longer stellar lifetimes causing a decrease in the super- 
nova rate. We note that in all models the wind power is more 
important than commonly recognized (e.g. in the study of su- 
pernova driven bubbles). This means that for stellar clusters, the 
supernova shells expand into a pre-exis ting cavity, rather th an in 
a dense star-forming ISM. For example fCho & Kangl d2008l) find 
that the volume of a supernova bubble is ~ 2 - 3 times larger and 
also hotter when the winds are taken into account, even though 
the wind power assumed in their study is very low compared to 
our estimates. 

The flow of wind and supernova ejecta inside cavities can 
be very complex (Mac Lo w et alj 120051) . and the propaga- 
tion is likely dominated by turbulent diffusion from magnetic 
field irregularities caused by the st ellar winds and supernovae 
(Balsara & Kiml2005tlParizotl 2004). Part of the gas will be ther- 
malized near the stellar association due to wind-wind collisions 
or by a termination shock against the turbulent medium inside 
the cavity , and this can be observed as a hot X -ray emitting 
plasma dTownslev et al J 120031: iGiidelet al.ll2008l) . However the 
majority of the mass will expand into the low density cavity. The 
effective propagation velocity is very uncertain, a nd is expected 
to be in the range 100 - 1000 km s -1 (for example Balsara et al. 
120081 finds velocities of -200 km s _1 ). 26 Al and 60 Fe can be im- 
portant for measuring this velocity, as their lifetimes are similar 
to the time it takes for the ejecta to cross a cavity blown by a 
young stellar cluster. If the velocity is high, these elements will 
quickly reach the wall of the cavity, and the brightness distri- 
bution will be given by the geometry of the walls. For slower 
propagation the fresh ejecta will be distributed inside the cavity, 
and the measured line widths will reflect the turbulence inside 
the cavity directly. 

5. Conclusions 

We have developed a populations synthesis code for the study 
of Galactic OB associations. The main aim of this work is the 
prediction of the output of energy, gas, 26 Al, 60 Fe and ionizing 
photons of a population of massive stars. We describe the 
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Fig. 13. The time profiles of the kinetic power, 26 Al, 60 Fe and 
UV, with the contributions from 4 mass ranges. The solid lines 
are for the full range of masses (same as the solid lines in Figs. 
SJandGI 



depende nce on the input physics , especially rotating stellar 
models dMeynet & Maederl l2005t iPalacios et all l2005h . We 
show that energy and matter output depend on the choice 
of stellar evolutionary and supernova models, whereas other 
parameters, such as the speeds of the winds and the stellar 
atmospheres only play secondary roles. This is a problem for the 
correct modelling of star-forming regions, and also means that 
comparison with observations (for example the y-ray emission 
from radioactive isotopes) can potentially yield important 
constraints for stellar and supernova models. Also the statistical 
variations for small numbers of stars in individual regions are 
addressed, which is important when interpreting observations of 
nearby, small star-forming regions. 

Our study shows that the ejected 26 Al and wind power from a 
population of massive stars is strongly enhanced by the effects 
of stellar rotation. Despite the recent downward revision of 
the mass-loss rates from massive stars by a factor of 2-3, the 
amount of 26 Al and the importance of the wind power is actually 
increased by including stellar rotation. 
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